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I^H . Abstract 

Time-dependent unitary transformations are used to study the Schrodinger equation for explicitly time- 
dependent Hamiltonians of the form H(t) — R(t) ■ J, where R is an arbitrary real vector-valued function of 

' time and J is the angular momentum operator. The solution of the Schrodinger equation for the most general 

Hamiltonian of this form is shown to be equivalent to the special case R = (1, 0, v{t)). This corresponds to the 
problem of driven two-level atom for the spin half representation of J. It is also shown that by requiring the 
magnitude of R to depend on its direction in a particular way, one can solve the Schrodinger equation exactly. 
y—\ ' In particular, it is shown that for every Hamiltonian of the form H(t) — R(t) ■ J there is another Hamiltonian 

with the same eigenstates for which the Schrodinger equation is exactly solved. The application of the results 
to the exact solution of the parallel transport equation and exact holonomy calculation for SU(2) principal 

■ bundles (Yang-Mills gauge theory) is also pointed out. 

Pk 
^— > ■ 

1 Introduction 
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In non-relativistic quantum mechanics the dynamics of pure states is determined by the Schrodinger equation, 

Hi/> = # , (1) 

where H is the Hamiltonian, ip is the state vector representing the state, the dot denotes a time-derivative, and 
Ti is set to unity. In general the Hamiltonian may be explicitly time-dependent, in which case the exact solution 
of the Schrodinger equation is in general not known. In terms of the time-evolution operator U = U(t) defined by 
U(t)ip(0) := ip(t), the Schrodinger equation (|l|) is written as 

H(t)U(t) = iU(t), 17(0) = 1. (2) 

An alternative expression for this equation is U(t) — Texp[— i f H(t')dt'], where T denotes the time-ordering 
operator. The purpose of this article is to derive some general sufficiency conditions to obtain the exact solution 
of Eq. (||) for the dipole Hamiltonians: 

3 

H(t)=Y,R a (t)Ja = R-J, (3) 

a=l 

where R a are real functions of time which do not simultaneously vanish, and J a are generators of the group SU(2) 
in some irreducible representation. 
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Following the same line of reasoning as in the Hamilton- Jacobi theory of classical mechanics, one can view 
the inverse ?7 _1 (t) = U'(t) of the evolution operator U(t) as a time-dependent quantum canonical transformation 
which sets the Hamiltonian to zero. In order to see this more clearly, let us first recall that in quantum mechanics 
the role of canonical transformations is played by unitary transformations of the Hilbert space. Now consider 
an arbitrary time-dependent unitary transformation ^(t) — ► i>'(t) = U(t)i])(t). Requiring this transformation to 
preserve the form of the Schrodinger equation ([!]) , one has 

H(t)^H'(t) = U{t)H{t)U\t)-iU{t)U\t) , (4) 
U(t)^U'(t) = U{t)U{t)U\Q) . (5) 

Hence, U'(t) induces a particular quantum canonical transformation which renders the transformed Hamiltonian 
zero. In other words, if one views the effect of a quantum canonical transformation as a change of frame in the 
Hilbert space, then the transformation induced by U^(t) takes one to a moving frame in which the state vector is 
stationary, i.e., ip'(t) = ip'(0). 

As it is manifestly seen from (^), quantum canonical transformations also resemble the non-Abelian gauge 
transformations of particle physics. Therefore, in a sense solving the Schrodinger equation (^|) is equivalent to 
finding an appropriate gauge in which the state vector is stationary.^] 

In this paper I shall try to demonstrate the utility of this simple observation in solving the Schrodinger equation 
for a large class of Hamiltonians of the form (^) . The basic idea pursued in this paper is to find a series of unitary 
(gauge) transformations which simplify the form of the Hamiltonian and yield previously unknown exactly solvable 
cases. Of course a priori there is no systematic method of choosing appropriate gauge transformations. However it 
turns out that at least for the systems considered here, one is guided by basic group theoretical properties of angular 
momentum operators and methods of quantum adiabatic approximation. Probably the most notable feature of 
this method which makes it so effective is its non-perturbative nature. 

2 Dipole Hamiltonian and Parallel Transportation in SU(2) Bundles 

Consider the Dipole Hamiltonian (^) 

H = H[R] =R J = r(sin0cos<M+sm0sin^2 + cos0J 3 ) = W(6, tp) J 3 W\6, tp) , (6) 

which describes the dynamics of a magnetic dipole in a changing magnetic field. Here R := (i? 1 , R 2 ,R 3 ) = (r, 8, tp) 
corresponds to the magnetic field vector expressed in units in which the Larmor frequency is set to unity, (r, 6, ip) 
are spherical coordinates, and 

W(6, tp) := e-^e-^V^ _ ( 7 ) 

Then an arbitrarily changing magnetic field corresponds to a curve C : [0, T] — > H 3 , R — R(t) — C(t). 

An application of the dipole Hamiltonian (|^) is in the parallel transportation in SU (2) principal fiber bundles 
(Yang-Mills theory). This is easily seen by recalling that parallel transportation |IJ] is defined in terms of a Lie 
algebra- valued one- form (gauge potential) A = A^Jadx^ according to 

g[C] =Ve^fc A =Ve- i tt? A ^ = Tg -< W)K<. ; (g) 

where C : [0, T] — > M is a curve in the base manifold M of the bundle (spacetime in Yang-Mills theory), and t is 
an arbitrarily chosen parameter of the curve C. It is very easy to recognize the last expression on the right hand 
side of (||) as the time-evolution operator U(T) for a Hamiltonian of the form (|^) with R a = x^(t)A a ^[x(t)]. Thus 
an exact solution of the Schrodinger equation for (|6|) yields as a special case the solution for the problem of parallel 
transportation and in particular the calculation of the holonomy elements and Wilson loop integrals in Yang-Mills 
theory. 

1 Note however that here there does not exist an analog of a non-Abelian gauge symmetry unless one deals with peculiar constraint 
systems such as those encountered in quantum cosmology. 
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Let us next recall the basic properties of the most general Hamiltonians of the form (|^) , |j . It is not difficult 
to see that the eigenvalues E n [R] and eigenvectors \n; R) of H[R] are given by: 

E n [R] = E n (r, 0, ip) = E n (r, 0, 0) = nr , with n = -j, -j + 1, • • • , j , (9) 
\n;R) = \n;(r,6,cp)) = \n;(r ,6,L P ))=W(6, ( p)\n;(r ,0,0)) , 6 € [0, tt), <p € [0, 2tt), (10) 

where j corresponds to the spin j-representation of SU(2) and determines the Hilbert space, and (tq, 0q, tpo) := 
(r(0),6(0),ip(0)). Hence, the Hamiltonian H is non-degenerate for r ^ 0. In order to avoid the complications 
caused by the sudden collapse of all the energy eigenvalues which occurs at r = 0, I shall only consider the case 
where the curve C does not pass through the origin, i.e., C(t) = R(t) em 3 - {0}. 

Note that {|n; (t, 0, </?))} forms a single- valued orthonormal basis of the Hilbert space for all ip € [0, 2tt) and 
€ [0, 7r) and that n; (r, 0, 0)) are the eigenvectors of H(r,0 — 0, ip = 0) = rJ 3 , i.e., J 3 |n; (r, 0, 0)) = n\n; (r, 0, 0)). 
For = t:, n; (f, 0, tp)) are not single- valued. This is due to the fact that the spectral bundle over H 3 — {0}, Qj, 
which yields |n; (t, 0, p)) as its local basis sections is not trivial. In the parameterization of M 3 — {0} used here the 
negative z-axis (0 = tt) is not included in the patch over which \n; (t, 0, ip)) are well-defined. To treat the negative 
z-axis, one must switch to new coordinates R 1 := (r' = r, 0' = tt — 0, p' = p). The eigenvectors |n; R') will then be 
single-valued everywhere except on the positive z-axis. In the following, I shall assume for simplicity but without 
loss of generality that the curve C does not intersect the negative z-axis. In the general case where C intersects 
the negative z-axis, one must make appropriate U(l) gauge transformations which relate \n;R) and \n;R'), 0. 

I shall also assume that |n; (0, 0)) and therefore \n; R) are eigenvectors of the total angular momentum operator, 
i.e., the Casimir operator \ J\ 2 = 5Z a=1 This is always possible unless R a are also quantized, |!|. The latter 
case will not be considered in the present article. 



3 Adiabatic Approximation and Reduction to Two-Dimensions 

In order to implement the idea of successive quantum canonical transformations, I shall begin using the results of 
the adiabatic approximation. One knows from the standard arguments of Born and Fock J^] and Kato || , that if the 
time-dependence of the Hamiltonian is adiabatic, then in time the eigenstates of the initial Hamiltonian H[R(0)] 
evolve into the eigenstates of the Hamiltonian H[R(t)]. This is actually very easy to see if one differentiates both 
sides of the eigenvalue equation 

H(t)\n;t) = E n (t)\n,t) , (11) 
and computes the inner product of both sides of the resulting equation with |m; t) for some m ^ n. This yields 

/ i d i , (m:t\H(t)\n\t) , . „. 

A mn := (m:t — \n:t) = x ' ' wl '/ , m^n. (12) 
V ' ] dt l ' ' E n (t)-E m {t) ^ V ' 

In Eqs. @ and @, H(t) := H[R(t)}, \n;t) := \n; R(t)), and E n (t) := E n [R{t)}. The adiabatic approximation is 
valid if and only if the right hand side of (|l^) is negligible. Now let us choose ip(0) — \n; 0), then in view of jl^), 
it is easy to show that ip(t) = e lotn ^' |n; t) does solve the schrodinger equation provided that 

a n (t) := 6 n (t) +7n(i) , 8 n (t) := - / E n (t')dt' 7n (t) := i f A nn {t')dt' . (13) 



The phase angles a n (T), S n (T), and 7«(T) for a closed curve C are known as the total, dynamical, and adiabatic 
geometrical (Berry) phase angles, 0. 

The adiabatic approximation which also includes the geometric phase effects, corresponds to approximating the 
time-evolution operator U(t) with 

U (t) :=Ye la "\n;t)(n;0\ . (14) 



In general the approximation U ~ Uq is not valid. However, one can compute Uq in terms of the eigenvalues and 
eigenvectors of the Hamiltonian and use Uq to perform a quantum canonical transformation. In the remainder of this 
section I shall show that indeed this canonical transformation simplifies the form of the Hamiltonian considerably. 
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In order to do this one must first calculate the matrix elements A mn which enter the calculation of a n and 
especially the term LiU^ in Eq. (Q) with U = Uq. This rather lengthy calculation leads to 

m(l - cos0)8 mn + - sin6» (e lip C m S mn -i + e~ lip C n 5 m - ln ) b + 

{e i<p C m 5 mn ^ - e-*C n 8 m -i n ) 9 , (15) 
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where C rn := y (j — m)(j + m + 1) = C_ m _i, and extensive use is made of the properties of J a and J± := J\ ±i J2, 
particularly 

e -i/3 J a Jb e tpj a _ cog ^ j b + £afcc sin /3 J c , a/t, 

J±|m;(r,0,0)> = a C ±m \m ± 1; (r, 0, 0)) , 

where e a (, c are components of the totally anti-symmetric Levi Civita symbol, with 6123 = 1. Furthermore, one can 
easily show that a n — na, S n — nS, j n = nj, where a = S + 7, and 

r(t') dt' , 7 = - / [1 - cos 9(t')} ip(t') dt' . (16) 
Jo 

These relations are then used to write down the expression for Uq, namely, 

U (t) = W(6(t),cp(t))e ia W J *Wi(6 ,<p ) , (17) 

where W is defined in Eq. (^). 

Next let us set Li — Uq in Eq. ([|) . Then using Eq. ( |l5| ) , one finds the expression for the transformed Hamiltonian 

Ho(t) = ^W(9 ,ipo)[Q(t) J+ + Q*(t) J-]Wi(6 ,(f>o) , (18) 

where 

fi(t) := e -i[«(*)+v(t)] [ S i n 6»(£) 0(f) + i0(t)] . (19) 

One can easily see that if 60 = <Po = 0, then W(9 ,(p ) — 1 and the expression (|l8| ) for the transformed 
Hamiltonian simplifies considerably. Hence, it is convenient to choose the coordinate system in such a way that 
$o = ^0 = 0, i.e., i?(0) = (0,0, ro), or alternatively make a further constant unitary transformation using Li = 
W* (9o,tpo) which leads to the Hamiltonian 

Hi(t) = uj(t) [coscr(i)Ji - sin a(t)J 2 ] , (20) 

where 51 =: we", i.e., 



uj(t) := \J9 2 + sin 2 9 b 2 , cr(t) := -a - ip + £ mod27r, 

sin 9 b . ^ 9 
cos£ := , sin( := — 

One can also combine the two unitary transformations by transforming H by Li = U\{t) witht/i(f) := Uo(t)W(9o, <po)- 

The Hamiltonian ( po| ) describes the dynamics of a magnetic dipole in a time-dependent magnetic field which is 
confined to the x-y plane, i.e., a Hamiltonian of the form @ corresponding to a planar curve Ci : [0,T] ->Er-{0}. 
Hence, the canonical transformation induced by U\ reduces the three-dimensional problem to a two-dimensional 
one. 



4 Exactly Solvable Cases 

Consider the Schrodingcr equation for the Hamiltonian H\. If the angular variable a happens to be constant, then 
this equation can be easily integrated. This is simply because in this case H± at different times commute and the 
transformed evolution operator is obtained by its exponentiation, i.e., 

rr//.\ _ -U(t)[cos a Ji -sin a ,J 2 ] 
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where £{t) := f*u(t')dt', 



.0 := .(0) = -<po + f(0) = + tan- ^ Q > m ] = + ^[-^] , (22) 
and & := d8/d<p. 

Having found the evolution operator U' for H\ , one can use Eq. (0) to write down the solution of the original 
Schrodinger equation (Q). This yields 

U(t) = U 1 (t)U' (t)Ul(Q) = U (t)W{9 ,ipo)U'{t)W^9 ,ip ), for a(t) = a . (23) 

Note that the parameters <7o and £ which enter the expression for U(t) are geometric quantities associated with 
the projection C of the curve C onto the unit sphere centered at the origin. In particular, I is the length of C. 
Furthermore for those portions of the curve C which project to a single point for an extended period of time, u> 
and consequently Hi vanish. This is reminiscent of the known fact that the adiabatic approximation is exact when 
the eigenvectors of the Hamiltonian are stationary. 

Another way of arriving at the same conclusion is by performing another quantum canonical transformation 
with U = lit := e~ lCT (') j3 . This leads to the transformed Hamiltonian 



H 2 =uj(t)Ji+&{t)J 3 . (24) 

Clearly for a = const, the Schrodinger equation for Hi is exactly solvable. Making a further canonical transforma- 
tion with U=Xl\:= e tt W Jl , one obtains 

H 3 = &[cos£(t) J 3 + sinl(t) J 2 ] . (25) 

which vanishes identically for a = const. Therefore, as expected the combined transformation U = (t/it^t^)^ leads 
to a frame in which the Hamiltonian vanishes and the state vector is stationary. Hence, the original time-evolution 
operator is given by U — U1U2U3. 

Let us next re-express the condition a = const, in terms of the original variables. Requiring a — 0, one finds 
the equivalent condition: r(t) = r*(t), where 



g c g ) 

/ ,\ n . dt v sin 9 (b * 

r*(t) := cosO (p 



Therefore, one has: 



dip V sin 9 > 

1 4. fX)2 
^ ^ v sin 6 I 



(26) 



Lemma 1: The exact solution of the Schrodinger equation ^) is given by provided that the magnitude 
of the magnetic field depends on its direction according to r{t) = r*(t). 

This is quite remarkable, for it indicates that for every Hamiltonian of the form @ for which r* does not vanish for 
extended periods of time, there exists another Hamiltonian with the same eigenvectors^] whose Schrodinger equation 
is exactly solvable. Note that for time intervals during which r» < 0, one can consider the time- reversed system 
where r* > 0. The evolution operator obtained for the time-reversed system yields the original time-evolution 
operator upon inversion. This leaves only the cases where r* vanishes, i.e., either tp = or 8' = sin#tan[sin# + c] 
for some constant c. A simple case where the latter equation is satisfied is 8 = ir/2 and c = — 1. This means that 
for the planar curves with 8 — it/ 2 such as C\, one cannot enforce the condition r = r* and the exact solution 
cannot be obtained in this way. Therefore a direct repetition of the same procedure for the Hamiltonian Hi will not 
lead to the exact solution. In the remainder of this section I shall demonstrate, however, that by a straightforward 
redefinition of the time one can generalize Lemma 1 further. 

Let us first note that for the case where u> = the exact solution is given by the adiabatic approximation. 
Hence, without loss of generality one can restrict to the case w^fl. In this case the length £ of the projection C of 

2 Note that the eigenvectors only depend on the direction of the magnetic field. 
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the curve C is a monotonically increasing function of time t. Therefore it can be used to parameterize the evolution 
of the system, i.e., replace t. Changing variables from t to I in the Schrodinger equation for the Hamiltonian H\ 
and making use of a; ^ 0, one has 

H l {l)Ui{l)=i± Ui{l), (27) 

where 

H X {1) cosa(£)Ji - sina(£)J 2 = e la{e)j3 J ie ^ ff(£)j3 . (28) 

This reduces the problem to the case of a magnetic field which traces a circular path in the x-y plane with an 
angular frequency, v :— da/d£ = (r — r„)/uj. Note that the presence of ui(t) on the right hand of Eq. ( p0| ) is quite 
essential in the redefinition of time. 

Let us next transform to the rotating frame defined by U = U^W := e - t(J W J 3 . I n view of Eq. (Q), this leads to 
the transformed Hamiltonian 

H 2 (£) = Jx+u(£)J 3 , (29) 

which describes a magnetic field with a constant rr-component and a variable z-component. Such systems are 
widely encountered in the study of nuclear and optical magnetic resonance. For a recent study of an iterative 
solution of the Schrodinger equation for this Hamiltonian see Refs. || [|. 

Note that for v = i>q = const., H2 is constant. Hence, the transformed time-evolution operator is given by 
{/3(f) := exp[— i£(Ji + V0J3)], and one has 

U(t) = f/i (t)02 (t))^ 3 (^(t))^ 2 + (0)^ (0) = U (t)W(9 , ip o )U 2 (£(t))U 3 (£(t))uU0)W^(9 , tpo) . (30) 

This concludes the derivation of the exact solution of the schrodinger equation for the case where a(£) = Co + vq£, 
alternatively, r(t) — r*(£) + vq uj(t). This is a generalization of Lemma 1. It states that even for the time periods 
during which = 0, the above procedure still leads to exactly solvable Schrodinger equations. More precisely, the 
following lemma holds. 

Lemma 2: The exact solution of the Schrodinger equation (^) is given by Eq. $3t\), provided that the 
magnitude of the magnetic field depends on its direction according to r(t) — r*(<) + VQLu(t), for some constant 
u . 

A direct consequence of this result is 

Corollary: For every Hamiltonian of the form there exists another Hamiltonian with the same eigen- 
vectors for which the Schrodinger equation is exactly solvable. 



5 Conclusion 

In this paper I have used a variety of time-dependent unitary transformations of the Hilbert space to obtain the 
exact solution of the schrodinger equation for a large class of explicitly time-dependent dipole Hamiltonians. This 
involved redefinition of the time variable which was a consequence of transforming to a moving frame via the 
inverse of the adiabatically approximate time-evolution operator. In this frame the natural choice for the evolution 
parameter turned out to be the length of the projection of the curve C traced by the tip of the magnetic field onto 
the unit sphere centered at the origin. 

The reduction of the general problem to that of the Hamiltonian H2 = J\ + v{£)J 3 may also be used to set 
up an approximation scheme for large u. This is due to the fact that v = (r — r„)/uj may be neglected for large 
u>, in which case Lemma 2 provides the solution. This is particularly effective for the dipole Hamiltonians which 
correspond to a planar curve C, for which r* = 0, e.g., Hi. For these Hamiltonians, the approximation is valid if 
the parameter r(t)/u>(t) is negligible. Note also that for such Hamiltonians if r(t) and ui(t) are proportional, then 
Lemma 2 yields the exact solution to the Schrodinger equation. 

Moreover, by successive application of the method used in this reduction, i.e., by replacing the original Hamilto- 
nian H by Hi and repeating the same analysis, one obtains an iterative solution of the Schrodinger equation which 
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yields a product expansion of the time-evolution operator. The condition of the termination of this expansion after 
a finite number of iterations may seem to lead to (possibly) more general exactly solvable cases. It turns out that 
this is in fact not the case. This is because enforcing the condition that the above expansion be terminated after 
the second iteration leads to u> = const. , which is certainly not more general than the conditions of Lemma 1 and 
Lemma 2. This marks a unique property of the Hamiltonians of type H(t) = J% + R 3 (t) J3. 

The results of this paper have direct applications in the computation of the holonomy elements and Wilson 
loop integrals in Yang-Mills theory where the gauge group is 5'C/(2).^| In this case the original parameters R a of 
the Hamiltonian (^|) are identified with i'M", where (A^) corresponds to the local connection one-form (gauge 
potential) and the gauge transformations correspond to quantum canonical transformations of the associated Hamil- 
tonian. Another area of application of the results of this article is in the calculation of non-Abelian U (2) geometric 
phases pj such as those encountered in the study of the three-level systems, Q |. 
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